Universal Thermoelectric Effect of Dirac Fermions in Graphene 
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We numerically study the thermoelectric transports of Dirac fermions in graphene in the presence 
of a strong magnetic field and disorder. We find that the thermoelectric transport coefficients 
demonstrate universal behavior depending on the ratio between the temperature and the width of 
the disorder-broadened Landau levels(LLs). The transverse thermoelectric conductivity a^y reaches 
a universal quantum value at the center of each LL in the high temperature regime, and it has 
a linear temperature dependence at low temperatures. The calculated Nernst signal has a peak 
at the central LL with heights of the order of fc_g/e, and changes sign near other LLs, while the 
thermopower has an opposite behavior, in good agreement with experimental data. The validity 
of the generalized Mott relation between the thermoelectric and electrical transport coefficients is 
verified in a wide range of temperatures. 
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Graphene has attracted enormous interest due to its 
unique electronic properties associated with the two- 
dimensional (2D) Dirac-fermion excitations [ij, includ- 
ing its thermoelectric properties 0^ • The thermopower 
and Nernst coefficient, measuring the magnitude of the 
longitudinal and transverse electric fields generated in re- 
sponse to an applied temperature gradient, are very sen- 
sitive to the semimetal nature of graphene. For instance, 
a large thermopower is expected near the Dirac point as 
Sxx ^ T/Ep while T and Ep are the temperature and 
the Fermi energy, respectively. This has been observed 
in recent experiments 0^3; where the maximum value 
of thermopower reaches dOfiV/K at T w 300i^. 

In a magnetic field, the electronic states of graphene 
arc quantized into Landau levels (LLs), as in the 2D 
semiconductor systems displaying the integer quantum 
Hall effect (IQHE) Q. However, the Hall conductiv- 
ity of graphene obeys an unconventional quantization 
rule axy — 4(n + l/2)e^/ft-, where n is an integer 
911. For the conventional 2D IQHE systems, theories [K 
13| predict that when the thermal activation dominates 
the broadening of LLs, all transport coefficients are uni- 
versal functions of Ep/fkOc and ksT/fiuic, where is 
the LL quantization energy. Sxx shows a series of peaks 
near the LL energies with height \n2{kB /e)/{n + 1/2), 
which is independent of the magnetic field or tempera- 
ture. The Nernst signal Sxy oscillates about zero near the 
LLs and enhances as the strength of the impurity scat- 
tering increases. The experimental results on graphene 
agree with these asymptotic behaviors except at the cen- 
tral LL, where Sxy has a peak instead with maximum 
value about 40iJ,V/K while Sxx becomes oscillatory 0- 
0] ■ The unusual behavior of Sxy and Sxx near the central 
LL has not been understood. 

While thermoelectric transports depend crucially on 



impurity scattering as well as thermal activation, the 
study of disorder effect on thermoelectric transports in 
graphene is still lacking. Another important question 
is to what extent the well-known Mott relation between 
thermoelectric and electrical transport coefficients [cf. 
Eq.([3])] is applicable for this system. In this Letter, we 
carry out a numerical study to address all the above is- 
sues. We show that thermoelectric transport coefficients 
are universal functions of the ratio between the temper- 
ature and the disorder-induced LL width, and display 
different asymptotic behaviors in different temperature 
regions, in agreement with experimental data. Our study 
also reveals that the distinct behaviors in the central LL 
is an intrisic property of the Dirac point where both par- 
ticle and hole LLs coexist. Furthermore, the generalized 
Mott relation is shown to be valid for a wide range of 
temperatures. 

We consider a rectangular sample of a 2D graphene 
sheet consisting of carbon atoms on a honeycomb lat- 
tice 1^. Besides the Anderson- type random disorder 
considered in Ref. 1J|, we also model charged impurities 
in substrate, randomly located in a plane at a distance 
d, either above or below the graphene sheet with a long- 
range Coulomb scattering potential. The latter type of 
disorder is known [l^ to give a more satisfactory inter- 
pretation of the transport properties of graphene in the 
absence of magnetic fields. When a magnetic field is ap- 
plied perpendicular to the graphene plane, the Hamilto- 
nian can be written in the tight-binding form 



{ij)a itr 
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where (cio-) creates (annihilates) a tt electron of spin 
a on lattice site i. t is the nearest-neighbor hopping in- 
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tegral with an additional phase factor aij due to the ap- 
phed magnetic field B. The magnetic flux pei^hexagon 
(j) — X^o "^u ~ W "^ith M an integer !l4]. For 
Anderson- type disorder, Wi is randomly distributed be- 
tween [—W/2, W/2] with W as the disorder strength. For 

charged impurities, Wi = -^J^a^/ V i^i - Ka)"^ + d"^ , 
where Ze is the charge carried by the impurities, k is 
the effective background lattice dielectric constant, and 
Ti and Ra are the planar positions of site i and impurity 
a, respectively. All the properties of the substrate (or 
vacuum in the case of suspended graphene) can be ab- 
sorbed into a dimensionless parameter — Ze'^ / (kTlvf), 
where vp is the Fermi velocity of the electrons. For sim- 
plicity, in the following calculation, we fix the value of 
distance d — I, randomly distribute impurities with the 
density as 1% of the total sites, and tune rg to control 
the impurity scattering strength. 

In the linear response regime, the charge current in 
response to an electric field or a temperature gradient 
can be written as J = ctE -|- «(— VT), where a and a are 
the electrical and thermoelectric conductivity tensors, re- 
spectively. These transport coefficients can be calculated 
with Kubo formula once we obtain all the eigenstates 
of the Hamiltonian (in our calculation, axx is obtained 



based on the calculation of the Thouless number [16|). In 
practice, we can first calculate the T = conductivities 
Gji^Ep), and then use the relation [Til 




FIG. 1: (color online). Zero temperature a^y and (Jxx as 
functions of the normalized Fermi energy near the Dirac point 
with Wl as the width of the central LL. (a) and (b) compare 
results for short-range random potential (with various W) 
and long-range Coulomb potential (with various r^), where 
the system size A*' = 96 x 48 and the magnetic flux = ||-. 
(c) and (d) compare results for three different strengths of 
magnetic flux, with the same system size A*" = 96 x 48 and the 
same impurity conflguration, where a long-range scattering 
potential with = 0.4 is assumed, (e) and (f) are for three 
different system sizes with = 0.3 and (f) — 
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aji{EF,T) = J deaji{e) 
aj^{EF,T) = j deaJ^{e){e- Ef) 
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to obtain the finite temperature electrical and ther- 
moelectric conductivity tensors. Here f{x) = 
'^l^^{x-EF)/kBT _|_ 2] is the Fermi distribution function. 
At low temperatures, the second equation can be ap- 
proximated as 



aji(EF,T) 



3e 



de 



(3) 



which is the generalized Mott relation [12, [l3 ■ 

In FiglU we show the calculated Hall conductivity a^y 
and longitudinal conductivity a^x at T = as functions 
of the Fermi energy near the Dirac point. From Fig. [T] 
(a) and (b), we observe that the results for the two dif- 
ferent types of disorder are very similar. The Hall con- 
ductivity exhibits two well-quantized plateaus —2e^/h 
and 2e^/h, being consistent with the quantization rule 
o'xy — 4(rt + I/2)e^//i. The direct transition between 
the two plateaus is accompanied by a pronounced peak 
in the longitudinal conductivity axx with the maximum 
value 2e^/h [l3l- Remarkably, once we scale the energy 
with the width of the central LL (Wl), which is de- 
termined by the full-width at half-maximum of the axx 



peak, all results fall into a single curve, though there 
are small deviations for axx at the peak tails. This is 
rather in accordance with the scaling theory on the quan- 
tum Hall hquid to insulator transition [l^, as Wl corre- 
sponds to an energy scale where the electron localization 
length (correlation length) is comparable to the system 
size. The scaling is further tested for different magnetic 
fields and different system sizes, as shown in Figs. [Ijc-d) 
and Figs. [2e-f), respectively. We conclude that, when 
Ep, Wl ^ fi-Wc, o'xx and axy are universal functions of a 
single parameter Ep/Wl ■ It is noteworthy that the uni- 
versal curves shown in Fig. 1 are for disorder strengths 
smaller than the critical values, which is relevant to the 
experiments With stronger disorder strength, an 

insulating state may appear at the Dirac point |14{ . 

In Fig. [21 we show the results of thermoelectric con- 
ductivities axy and axx at finite temperatures. As seen 
from Figl^Ka) and (b) , axy displays a series of peaks near 
the LL energies, while axx oscillates and changes sign at 
the LL energies. These behaviors are similar to those in 
the conventional IQHE systems I4] , but some important 
differences exist. First, at low temperatures, the peak of 
axy at the central LL is higher and narrower than at the 
other LLs, which indicates that the impurity scattering 
has a different effect on the central LL and the rest LLs. 
Second, axy and axx are symmetric and antisymmetric 
about the Dirac point (zero energy) , respectively, as they 
are even and odd functions of the Fermi energy, rather 
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FIG. 2: (color online). Thermoelectric conductivities at finite 
temperatures, (a) and (b) show axy{EF,T) and axx{EF,T) 
as functions of the Fermi energy at different temperatures, 
(c) shows the temperature dependence of axy{EF,T) and 
oixx{Ef,T) for certain fixed Fermi energies, (d) compares 
the results from numerical calculations and from the gen- 
eralized Mott relation at two characteristic temperatures, 
UbT/Wl = 0.05 and UbT/Wl = 1. The other parameters 
are taken to be iV = 96 x 48, = 27r/48, and = 0.3 with 
Wl/I = 0.0195. 



tion still predicts the correct asymptotic behavior. This 
is due to the fact that the conductivities display uni- 
versal dependence on Ep/kBT, either a power law or 
an exponential function in different regimes, which can 
all be captured by the Mott relation. This proves that 
the generalized Mott relation is asymptotically valid in 
Landau-quantized systems, as suggested in Ref. [l^ . 
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FIG. 3: (color online). Universal behaviors of Oxx and axy 
near the central LL at fcsT/Wi = 1. (a)(b) show them as 
functions of the renormalized Fermi energy for three disorder 
strengths, and (c)(d) compare results for different magnetic 
fields. The parameters chosen are shown in the figure. 



than periodic functions. We also find that, depending on 
the relative strength between ksT and Wl, the thermo- 
electric conductivities show different universal behaviors. 
When ksT < Wl and Ep < Wl, we find that both 
axy and axx are linear in T (verified by a log-log plot 
to the lowest temperature accessible by numerics). This 
indicates that within the mobility edge where extended 
states dominate, the diffusive transports play important 
roles as in semiclassical Drude-Zener regime. When the 
Fermi energy falls deep inside the mobility gap or fc^T 
becomes comparable to or greater than Wl , thermal acti- 
vation dominates. In this regime, axx assumes the Arrhe- 
nius form (1/T)e^^^/'''^"^. Meanwhile, the heights of the 
peaks in axy for all LLs saturate to a value 2.77kBe/h, 
as seen from FiglJJc). This value matches exactly the 
universal number {\n2)kBe/h predicted for the conven- 
tional IQHE systems in the case where thermal activation 
dominates 12|, [l3| , with an additional degeneracy factor 



To examine the validity of the generalized Mott re- 
lation, we compare the above results with those calcu- 
lated from Eq.([3|), as shown in FigHJd). The Mott re- 
lation, which was historically derived from the semiclas- 
sical Boltzmann equation, is a low-temperature approx- 
imation and predicts that thermoelectric conductivities 
are linear in temperature. This is in agreement with our 
low-temperature results. At high temperatures, thermo- 
electric conductivities deviate from the linear-T depen- 
dence. However, if we take into account the finite tem- 
perature values of electrical conductivities, the Mott rela- 



Given the single-parameter scaling behaviors of the 
T — conductivities, it is straightforward to show with 
Eq. ^ that aji and aji at finite temperatures are univer- 
sal functions of WL/ksT and Ep/kBT, or 

CiEF,T) = CoSc{EF/WL,WL/kBT), (4) 

where C stands for one of the transport coefficients and 
Scix,y) is a universal function. We can directly verify 
the universal relations Eq.(|3|). Here, we pick a tempera- 
ture kBT ~ Wl, corresponding to 5'£(a:;, 1), and compare 
the results around the central LL for different disorder 
strengths and different magnetic fields, as shown in Fig. 
[3l Indeed, once we scale the Fermi energy with Wl, axx 
and axy as functions of Ep/Wl collapse into the same 
corresponding curves. 

We note that the asymptotic behaviors and the values 
of the peak heights of our calculated axx and axy are in 
agreement with the experimental results |3] . We further 
calculate the thermopower Sxx and the Nernst signal Sxy 



using [12 1 



ik^kj 7 



(5) 



k=x,y 



which are directly determined in experiments by mea- 
suring the responsive electric fields. The results for three 
central LLs are shown in Fig. |31 Sxy {Sxx) has a peak at 
the central LL (the other LLs), and changes sign near the 
other LLs (the central LL). The height of the Sxx peak at 
n = -1 LL is found to be 37^J.V/K for /cbT = 0.2Wl and 
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FIG. 4: (color online). Calculated (a) Sxx and (b) Sxy as 
functions of the Fermi energy for different temperatures. The 
parameters chosen here are the same as in Fig[2l 



56iJ,V/K for fcsT — 0.5Wl, which is in agreement with 
the maximum measured value 48/xT^/-ftr0]. This is also 
in agreement with the theory predication that, in the 
absence of disorder and at low temperatures, the peak 
value of Sxx is dominated by axy/oxy and takes a uni- 
versal value (2/3)ln2fcB/e « AQfiV/K. In the presence 
of disorder and at finite temperatures, the peak value is 
slightly bigger and the peak position is shifted toward 
Ep = 0. At Ep = 0, both a^y and a^x vanish, leading to 
a vanishing S^x- Around the zero energy, because (JxxOtxx 
and (JxyOixy have opposite signs, depending on their rel- 
ative magnitudes, Sxx increases or decreases when the 
Fermi energy is increased passing the Dirac point. In 
our calculations, we find that the second term is always 
dominant, which is different from the experimental ob- 
servation. This might be related to the unexpected large 
value of Gxx ~ 6e^//i observed in experiments. Figure |3] 
actually shows the result when we rescale axx to the ex- 
perimental value instead of the theoretical result 2e^//i. 
We can then obtain the same asymptotic behavior as in 
experiments for fc^T = Q.2Wl. On the other hand, S^y 
has a peak structure, which is dominated by a^yla^x- 
With the rescaled a^^^ we find that the peak height is 
40^F/-fC at ksT — 0.2Wl, which is comparable with the 
experimental value 20 — 40iJ,V/K. The distinct thermo- 
electric behaviors near the central LL can be traced down 
to the Berry phase anomaly of the Dirac point. The cen- 
tral LL in graphene in fact consists of two degenerate LLs, 
one from particles and another from holes, which are pro- 
tected by the particle-hole symmetry. While this makes 
no distinction in electrical transports when the energy 
is tuned through the LL, electrons and holes contribute 
differently in thermoelectric transports than in electri- 
cal transports. Here, the contributions of equal numbers 
of electrons and holes moving along the same thermal 
gradient direction cancel with each other in Sxx and are 
additive in Sxy Other LLs, without this property, have 
similar behaviors as in conventional IQHE systems. 



In summary, we have investigated the thermoelectric 
transports in graphene by a numerical study on the lat- 
tice model in the presence of both disorder and a mag- 
netic field and obtain results in agreement with experi- 
ments. 
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